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Abstract — There are several cases in wireless communications 
tlieory wliere the statistics of the sum of independent or corre- 
lated Nakagami-m random variables (RVs) is necessary to be 
known. However, a closed-form solution to the distribution of 
this sum does not exist when the number of constituent RVs 
exceeds two, even for the special case of Rayleigh fading. In 
this paper, we present an efficient closed-form approximation for 
the distribution of the sum of arbitrary correlated Nakagami- 
m envelopes with identical and integer fading parameters. The 
distribution becomes exact for maximal correlation, while the 
tightness of the proposed approximation is validated statistically 
by using the Chi-square and the Kolmogorov-Smirnov goodness- 
of-flt tests. As an application, the approximation is used to study 
the performance of equal-gain combining (EGC) systems oper- 
ating over arbitrary correlated Nakagami-m fading channels, 
by utilizing the available analytical results for the error-rate 
performance of an equivalent maximal-ratio combining (MRC) 
system. 

Index Terms — Nakagami-m fading, arbitrary correlation, ap- 
proximative statistics, equal gain combining (EGC), maximal 
ratio combining (MRC) 



I. Introduction 

THE analytical determination of the the probability distri- 
bution functions (PDF) and the cumulative distribution 
functions (CDF) of the sums of independent and correlated 
signals' envelopes is rather cumbersome, yielding difficulties 
in the theoretical performance analysis of some wireless 
communications systems (1). A closed-form solution for the 
PDF and the CDF of the sum of Rayleigh random variables 
(RVs) has not been presented for more then 90 years, except 
when the number of RVs equals two. The famous Beaulieu 
series for computing PDF of a sum of independent RVs were 
proposed in |2j. Later, a finite range multifold integral for 
PDF of the sum of independent and identically distributed 
(i.i.d.) Nakagami-ni RVs was proposed in |3|. A closed-form 
formula for the PDF of the sum of two i.i.d. Nakagami-?7i 
RVs was given in |4'l-|'6l. Exact infinite series representations 
for the sum of three and four i.i.d. Nakagami-yn RVs was 
presented in Q, although their usefulness is overshadowed 
by their computational complexity. 
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The most famous application, where these sums appear, 
deals with the analytical performance evaluation of equal 
gain combining (EGC) systems |8|-|13|. Only few papers 
address the performance of EGC receivers in correlated fading 
with arbitrary-order diversity. In |14|, EGC was studied by 
approximating the moment generating function (MGF) of the 
output SNR, where the moments are determined exactly only 
for exponentially correlated Nakagami-rn channels in terms 
of multi-fold infinite series. A completely novel approach 
for performance analysis of diversity combiners in equally 
correlated fading channels was proposed in |15|, where the 
equally correlated Rayleigh fading channels are transformed 
into a set of conditionally independent Rician RVs. Based on 
this technique, the authors in |16| derived the moments of the 
EGC output signal-to-noise ratio (SNR) in equally correlated 
Nakagami-m channels in terms of the Appell hypergeometric 
function, and then used them to evaluate the EGC performance 
metrics, such as the outage probability and the error proba- 
bility (using Gaussian quadrature with weights and abscissas 
computed by solving sets of nonlinear equations). 

All of the above works yield to results that are not expressed 
in closed form due to the inherent intricacy of the exact sum 
statistics. This intricacy can be circumvented by searching 
for suitable highly accurate approximations for the PDF of 
a sum of arbitrary number of Nakagami-m RVs. Various 
simple and accurate approximations to the PDF of sum of 
independent Rayleigh, Rice and Nakagami-?7i RVs had been 
proposed in llT7ll - ll2n . which had been used for analytical 
EGC performance evaluation. Based on the ideas given in fT|, 
the works lfT8l - ll2n use various alternatives of the moment 
matching method to arrive at the required approximation. 

In this paper, we present a highly accurate closed-form 
approximation for the PDF of the sum of non-identical ar- 
bitrarily correlated Nakagami-m RVs with identical (integer) 
fading parameters. By applying this approximation, we eval- 
uate the performance of EGC systems in terms of the known 
performance of an equivalent maximal ratio combining (MRC) 
system 1221 . Il24l . thus avoiding many complex numerical 
evaluations inherent for the methods presented in the afore- 
mentioned previous works for the EGC performance analysis. 
Although approximate, the offered closed-form expressions 
allow to gain insight into system performance by considering, 
for example, operation in the low or high SNR region. 



II. An accurate approximation to the sum of 

ARBITRARY CORRELATED NAKAGAMI-m ENVELOPES 

Let Z he a sum of L non-identical correlated Nakagami-m 
envelopes, {Zk}J^^i, defined as 



z = j2Zk 



(1) 



fc=i 



The envelopes {Zk}j^^i are distributed according to the 
Nakagami-m distribution, whose PDF is given by yj 



M (z) 



n. 
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m 



^ ^2 
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with arbitrary average powers E[Zl] — ilk, ^ 1^ k < L, and 
the same (integer) fading parameter m^. The power correlation 
coefficient between any given pair of envelopes (Zi,Zj) is 
defined as 



Prj 
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where E[-], cov(-, •) and var(-) denote expectation, covariance 
and variance, respectively. 

We propose the unknown PDF of Z be approximated by 
the PDF of R defined as 



R = 



\ 



k=l 



(4) 



where Rk, 1 < k < L, denote a set of L correlated 
but identically distributed Nakagami-m envelopes with same 
average powers, E[Rl] — flji, and same fading parameters, 
niR. The power correlation coefficients between any given pair 
{Ri, Rj) is assumed equal to that of the respective pair of the 
original envelopes {Zi,Zj), pij. 

The statistics of R^ is easily seen to be equal to the statistics 
of the sum of correlated Gamma RVs. Thus, the MGF of R^ 
is represented by ll26l Eq. (11)] 



Mji2 (s) = dct I 
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where I is the L x L identity matrix and A is the L x L 
positive definite matrix (denoted as the correlation matrix) 
whose elements are the square roots of the power correlation 
coefficients. 
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The L eigenvalues of the correlation matrix A are denoted by 
Afe, 1 < fc < L. 

Throughout literature, the PDF of R^ is determined by using 
several different approaches that result in alternative closed- 
form solutions, two of which are given by 1231 Eq. (29)] and 



ll24l Eq. (10)]. After a simple RV transformation, these two 
alternatives for the PDF of R are expressed as 
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where $2 (') is the confluent Lauricella hypergeometric 
function of L variables, defined in ||33l and ll24l Eqs. (9)- 
(10)]. Note that dHJ is here presented to demonstrate existence 
of an exact closed-form solution, whereas ^ is much more 
convenient for accurate and efficient numerical integration. For 
example, the PDF may be obtained using the Gauss-Legendre 
quadrature rule 01 Eq. (25.4.29)] over ^ ll23l . 

Next, we apply the moment matching method to determine 
the parameters flR and niR of the proposed approximation 
dT)-® to the PDF of Z. In wireless communications, moment 
matching methods are most typically applied to approximate 
distributions of the sum of log-normal RVs |29]. Most recently, 
a variant of moment matching, matching of the normalized 
first and second moments, had been applied to arrive at an 
improved approximation to the sum of independent Nakagami- 
m RVs via the a-p distribution ||20]| - ||2T1 . 

We arrive at required approximation by matching the first 
and the second moments of the powers of Z and R, i.e., the 
second and fourth moments of the envelopes Z and R, 

E[Z^] = E[R% E[Z'^] = E[R'^] . (9) 

Matching the first and the second moments of the powers aids 
the analytical tractability of the proposed approximation due 
to the availability of the MGF of R^ in closed form, given by 
(|5]l. The second and the fourth moments of R are determined 
straightforwardly by applying the moment theorem over (|5]l, 
yielding 
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Introducing dTOl l and ( fTTT i into (|9]l, one obtains the unknown 
parameters for the statistics of R 
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Using the multinomial theorem and IT] Eq. (137)], the second 
and the fourth moments of Z are determined as 

L 

= J2^k 



E[Z^ 



2r^{m, + 1/2) 



fc=i 

L 



i=l j=i+l 



L 

Y, ^/fm]2Fi{-l/2,-l/2;m,;p,j), (13) 






X2F1 (-l,-l;TO2;pij) + 



1=1 j=i+l 

4r(m2 + 3/2)r(m^ + l/2) 



r7?,^r2(?Ti2) 



2=1 J=Z+1 

L L L 

+ 12 E E E ^nv/fW^^I^m^^^j] 
m—1 i—m+1 j = z+l 
L L L 

+ 12 E E E \/^^^^V^E[Z^ZfZJ] 
rn—1 i—m-\-l j—i+1 
L L L 

+12 E E E \/^lJhn,E[Z^Z,Z]] 

m—1 i—m+1 j—i+1 



L L L L 

24 E E E E V^^^AJh^E[Z.mZ„Z,Z_ 

?n=l n—rn+1 i—n-\-l j—i-\-l 



3\ ' 



(14) 

where 2^1 (■) is the Gauss hypergeometric function ISTl . The 
joint moments £'[Z^Z,Zj], £'[^^^2^^], £'[Z„Z,Z|] and 
-E[ZmZ„ZiZj] are not known in closed-form for arbitrary 
branch correlation. Exact closed-form expressions are avail- 
able only for some particular correlation models, such as the 
exponential and the equal correlation models. For the case or 
arbitrary correlation, we utilize the method presented in fZfl, 
where an arbitrary correlation matrix A is approximated by its 
respective Green's matrix, followed by the application of the 
available joint moments of the exponential correlation model. 



A. Equal correlation model 

Equal correlation typically corresponds to the scenario of 
multichannel reception from closely spaced diversity antennas 
(e.g., three antennas placed on an equilateral triangle). This 
model may be employed as a worst case correlation scenario, 
since the impact of correlation on system performance for 
other correlation models typically will be less severe ll22l . 

mni. 

For this correlation model, the power correlation coefficients 
are all equal. 



Ptj^ p, i^ j, < p < 1 . 



(15) 



When TTiz is assumed to be integer, the unknown joint mo- 
ments in (O can be expressed in closed-form as |[T6l Eq. 

(43)] 
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where the coefficients W{- ■ ■) are determined as 
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with i^^ ' (•) denoting the Lauricella Fa hypergeometric func- 
tion of N vai-iables, defined by EH Eq. (9.19)] and El Eqs. 
(11)-(13)]. 

Note that the coefficient W{2,1,1) needs to be evaluated 
when L > 3, whereas the coefficient T4^(l, 1, 1, 1) needs to be 
evaluated when L > 4. In Appendix A, T4^(2, 1, 1) is reduced 
to the more familiar hypergeometric functions, attaining the 
form given by jA.lh W(l, 1, 1, 1) requires numerical evalua- 
tion of the Lauricella Fa function of 4 variables, which can 
be computed with desired accuracy by using one of the two 
numerical methods presented in 1251 Section IV. A]. 

The assumption of equal average powers, fi/j — ilz, 1 < 
k < L, yields independence of mji from flz- For this case. 
Table I gives the values of niji for several combinations of p, 
L and mz ■ The use of Table I aids the practical applicability 
of our approach for the case of equal average powers. 

For the equal correlation model, the eigenvalues of A are 
exactly found as Ai = 1 + (L — l)y/p and Afc = (1 — ^/p) for 
2 < fc < L, so the statistics of R^ is identical to that of the 
sum of a pair of independent Nakagami RVs. Thus, the MGF 
of i?2 is given by ESI Eq. (9.213)], whereas the PDF of R is 
given by EQI Eq. (9.208)] 
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where i-F'i() is the Kummer confluent hypergeometric func- 
tion EI] Eq. (9.210)]. 

B. Exponential correlation model 

Exponential correlation typically corresponds to the sce- 
nario of multichannel reception from equispaced diversity an- 
tennas in which the correlation between the pairs of combined 
signals decays as the spacing between the antennas increases 
L22J, Q0|. 

For this correlation model, the power correlation coefficients 
are determined as 



Pij = p' 



-j\ 



0< p<l 



(20) 



The unknown joint moments in ( fT4b . E[Z^ZiZj], 
E[Z^ZJZj], E[Z„,Z,Z'j] and E[Z^Z^Z,Zj], can be 
calculated from lfT4l Eqs. (11) and (12)]. The Appendix B 
derives simpler alternatives to lfT4l Eqs. (11) and (12)], which 



TABLE I 
Fading parameter m^ for some feasible scenarios with equal correlation 
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1.93 


1.9217 


2.943 


2.9258 


2.9168 


0.2 


0.9195 


0.8884 
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1.9102 
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1.8569 


2.9068 


2.8715 


2.8518 


0.4 


0.9156 


0.8841 
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1.907 


1.8722 


1.8535 


2.9039 
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0.9304 
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0.8929 


1.9242 


1.8971 


1.8831 


2.9222 


2.8944 


2.8799 
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0.9587 


0.9445 


0.9374 


1.956 
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1.9333 


2.9553 


2.9399 
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TABLE II 

Fading parameter m^ for some feasible scenarios with exponential correlation 
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involve a single infinite sum and a familiar hypergeometric 
function. 
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In (ISlT l. (?T.i,ri2,ri3) = (2,1,1) for the calculation of 
E[Z'^^ZiZj], (^1,^2,^3) = (1,2,1) for the calculation of 
E[ZmZfZj] and (711,71.2,71.3) = (1,1,2) for the calculation 
of E[Z„iZiZj]. The matrix A — [Sij] is the inverse of A's 
principal submatrix composed of the r7i-th, z-th and j-th rows 
and columns of A, whereas the matrix 'if = [ipij] is the 
inverse of A's principal submatrix composed of the 771-th, 77-th, 
7-th and j-th rows and columns of A, 

1 / / "■ ~^ 

J- Y Pni i v/ Prnj 

'Pirn ^ yPiJ 



The exactness of (l2ni-(l22b arise from the fact that both 
matrices A and ^ are tridiagonal matrices due to ( |20] i 1271 
Section IV]. Introducing (I2ni-(|22]| into ( fT4] i. one obtains the 
closed-form expression for E[Z*], which is omitted here 
for brevity. Combining (fT3]l-(fT4]i into ( fT2] i. one obtains the 
unknown parameters flu and niR for the statistics of R. 

The assumption of equal average powers, flk ~ ^z, 

1 < k < L, again renders independence of rriB. from flz for 

the exponential correlation model. Under such assumptions. 

Table II displays the values of mj^ for several illustrative 

(21) combinations of p, L and niz- 



C. Arbitrary correlation model 

In the general case of arbitrary branch correlations, the cor- 
relation matrix A is approximated by its appropriate Green's 
matrix, C, utilizing the method presented in 1271 Section 
IV]. Since principal submatrices of Green's matrices are also 
Green's matrices, the matrices A and ^, defined by ( i23l l. are 
determined to be tridiagonal, yielding direct applicability of 
the results presented in Section II. B to determine the unknown 
parameters iln and niji for the statistics of R. Thus, the 
statistics of Z are approximated by the statistics of R, whose 
arbitrary correlation matrix A is approximated by the Green's 
matrix C. In the following subsection, we illustrate the highly 
accurate approximation to the PDF of Z facilitated by this 
approach. 
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D. Validation via statistical goodness-of-fit tests 

We now statistically validate the proposed PDF approxima- 
tions for equal, exponential and arbitrary branch correlation by 
using two different goodness-of-fit tests. The Chi-square (C- 
S) and Kolmogorov-Smirnov (K-S) tests provide two different 
statistical metrics, Xn ^i^^ ^n, which describe the discrepancy 
between the observed samples of Z and the samples expected 
under the analytical distribution (JT])-©. 

Each metric is averaged over 100 statistical samples, where 
each statistical sample comprises of 10000 independent ran- 
dom samples of Z. The random samples of Z are generated 



TABLE III 
Significance levels of C-S and K-S tests for goodness 
OF fit between the exact and the approximative distributions of Fig. 1 
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0.02 


0.17 
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< 0.001 


< 0.001 


< 0.001 


< 0.001 
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0.04 


0.03 


0.2 


0.18 


< 0.001 


< 0.001 


< 0.001 


< 0.001 



TABLE IV 
Significance levels of C-S and K-S tests for goodness 
of fit between the exact and the approximative distributions of fig. 2 
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Fig. 1. Exact (obtained by simulation) and the approximative analytical 
PDFs to the sum of equally correlated Nakagami-m RVs with equal average 
powers, when Qz = 1 



Fig. 2. Exact (obtained by simulation) and the approximative analytical PDFs 
to the sum of exponentially con'clated Nakagami-m RVs with equal average 
powers, when Slz = 1 



by computer simulations of correlated Nakagami-m RVs based 
on the method proposed in 1281 Section VII]. 

For each metric, we calculate the significance level a from 
the C-S and K-S distributions, respectively denoted as acs 
and aKS- The significance level a represents the probability of 
rejecting the tested null hypothesis {Hq'. "the random samples 
of Z, obtained from ([T]), belong to the distribution given by (|7]i- 
dlll"), when it is actually true. The small values of a indicate 
a good fit. 

Note that, significance levels a less then 0.2 still indicate a 
good fit, due to the rigourousness of both C-S and K-S tests 
in accepting the null hypothesis Hq. 

1) Equal and exponential correlation: For the equal and 
exponential correlation models, the goodness-of-fit testing is 
conducted for combinations of the followings input parame- 
ters: L — 2 and 5, mz = 1 and 3, p = 0.2 and 0.7, whereas the 
average powers of Zk are assumed equal to unity (ftz ~ 1). 
The needed fading parameter m/j of distribution (I7]l-(l8]l is 
obtained directly from Tables I and II, whereas the average 
power ilji is calculated from (fT2T i. 



Figs. 1 and 2 depict the excellent (visual) match between 
the histogram obtained from generated samples of Z and the 
proposed approximation, for the cases of equal and exponential 
correlation models, respectively. Tables III and IV complement 
Figs. 1 and 2, by presenting the significance levels a for 
the corresponding input parameters' combinations. The Table 
III and the Table IV entries reveal the very low significance 
levels a for all input parameters' combinations, thus proving 
an excellent goodness of fit in statistical sense. 

2) Arbitrary correlation: For illustrative purposes, we use 
same two example correlation matrices from 1271 Sections V.B 
and VD], Sajin and I]4_circ, here denoted as Ai and A2, 
respectively. They are approximated by their Green's matrices 
C'ajin and C4_circ, here denoted as Ci and C2, respectively. 

Using Ci and C2, one obtains the needed tridiagonal 
matrices A and >& from their definitions given by (l23l l. The 
required joint moments are then calculated from ( I^Tl i and ( l22l i. 
which are then substituted into JTJI i and (O to calculate 
E[Z^] and E[Z% and then ([T2| is used to describe the 
statistics of R. 



TABLE V 

Significance levels of C-S and K-S tests for goodness of fit 

between exact and approximative distributions of fig. 3 
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m,z 


= 3 




acs 


OiKS 


acs 


OiKS 


Ai 


0.19 


0.18 


< 0.001 


< 0.001 


A2 


0.12 


0.11 


< 0.001 


< 0.001 



-A m =3 (approx) 
A^, m^ = 1 (approx) 

- A^, m^ = 3 (approx) 
Simulation (exact) 




Fig. 3. Exact (obtained by simulation) and the approximative analytical PDFs 
to the sum of correlated Nakagami-m RVs with equal average powers, when 
Qz = 1). 2nd congelation matrices Ai and A2 



from ( fTSl i. as 

L L 

f^fl = - ^ ^ y/ili^j, rriR 



{E[Z' 



i=i j=i 



E[Zi] - {E[Z'^]Y 



(25) 

respectively, where the latter equality is attributed to the 
definition of the Nakagami-TTi fading parameter, given by [T] 
Eq. (4)]. 

Thus, maximal correlation yields (l7]l-([8]l as an accurate 
distribution of Z, when our moment matching approach is 
appUed. This conclusion further validates our approach. 

III. Application to the performance analysis of 

EGC RECEIVERS 

We now consider a typical L-branch EGC diversity receiver 
exposed to slow and flat Nakagami-m fading. The envelopes of 
the branch signals Zk are non-identical correlated Nakagami- 
m random processes with PDFs given by (|2]l, whereas their 
respective phases are i.i.d. uniform random processes. Each 
branch is also corrupted by additive white Gaussian noise 
(AWGN) with power spectral density No/2, which is added 
to the useful branch signal. In the EGC receiver, the random 
phases of the branch signals are compensated (co-phased), 
equally weighted and then summed together to produce the 
decision variable. The envelope of the composite useful sig- 
nal, denoted by Z, is given by ([T]i, whereas the composite 
noise power is given by (t'^egc ~ LNq/2, resulting in the 
instantaneous output SNR given by 



Fig. 3 depicts the excellent (visual) match between the 
histogram obtained from generated samples of Z and the 
proposed approximation Q-®, for the two example corre- 
lation matrices Ai and A2. Table V complements Figs. 3, 
by revealing the very low significance levels a^, thus again 
proving an excellent goodness of fit. 



E. Validation in case of maximal correlation 

We now consider the case of maximal correlation coefficient 
between any pair of Nakagami-rn envelopes Zi and Zj, i.e., 
Pij = 1. It indicates a perfect linear relationship between these 
pairs, which, after applying the model from |[T6l Eq. (37)], can 
be defined as Zi — ^/TTi Zq for I < i < L, where Zq is an 
auxiliary Nakagami-m RV with unity average power and same 
fading parameter ttlz. After replacing the latter expression into 
([T]i, Z is transformed into a Nakagami-m RV with fading 
parameter niz and average power of 



1 



E[Z^] = 




L L 
i=l j = l 



(24) 



which agrees with (fTST l when pij -^ 1. 

Replacing pij = 1 into (|6]l, the L — 1 eigenvalues of the 
matrix A turn up equal to 0, except Ai — L. After plugging 
these eigenvalues into (|5]l, R is transformed into a Nakagami- 
m RV with fading parameter mji and average power Lftji. 
After the moment matching, iljj and niR can be obtained 



lEGC 



^" EGC 



LNo 



E^O - E^^ 



(26) 



\k=l 



\k=l 



where RVs Gk — Zk/^/LNo, 1 < k < L, form a set of L non- 
identical equally correlated Nakagami-m RVs with i?[G'|] = 
7fc/i, same fading parameters mz and correlation coefficient 
Pij between branch pair {i,j). Note that 7^ = Uk/No denotes 
the average SNR in fc-th branch. 

Using the results from Section II, the MGF and the PDF of 
( l26b can be approximated using (|5]l and ©-(O, respectively, 
when ilj^ is replaced by 7/j = ilj^/{LNo). These approxima- 
tions are then used to determine the outage probability -F-yEcc 
and the error probability Pegc of an L-branch EGC systems 
in correlated Nakagami-m fading with high accuracy. 

A. Outage probability 

The outage probability of the EGC system with arbitrary 
correlated Nakagami-m fading branches, whose output SNR 
drops below threshold i, is approximated by the known outage 
probability expressions of an equivalent MRC system 1231 Eq. 
(28)], m Eq. (13)], 



'^KEfe;oarctan(x%^) 
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ca sin 



xt 



dx 
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L-l 
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Equal correlation (approx) 

— Exponential correlation (approx) 
o Simulation (exact) 
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Equal correlation (approx) 

— Exponential correlation (approx) 
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Fig. 4. Exact and approximate en'or probabilities of an EGC receiver with Fig. 5. Exact and approximate error performance of an EGC receiver with 
correlated Nakagami-m branches, when m^ = 2, fi = and p = 0.7 correlated Nakagami-m branches, when m^ = 2, /i = 0.3 and p = 0.7 



r{l + Lmn) \nR 



rriR 



Lrriji 



1 



x$ 



(L) 



mR,mR,. 



det(A) 
,mjf, 1 + LniB-, 



rriR 



t. 



mR 



t. 



mR 



(27) 



For the equal correlation model, dZTJ l can be simplified using 
m Eq. (2.1.3(1))]. 

B. Average error probability 

The average error probability of the correlated Nakagami-m 
EGC system with BPSK modulation / coherent demodulation 
is approximated using the available expressions for the average 
error probability of the equivalent MRC systems. Based on 
Il30l Eq. (9.11)] and gH Eq. (17)], the error performance of 
this EGC system is alternatively approximated as 

1 r/2 

Pegc-bpsk ~ Pmrc-bpsk — — / Mr2 



-1 



de 

(28) 



T{LmR + 1/2) 
2^T{LmR + 1] 



mR 

nR 



Lin J 



1 



(imR 



det(A) 

+ l/2,mR,...,mR;LmR- 
m,R ruR mR 



^rM ^rM 



^R^L 



(29) 



In (|28ll, Mr2 (•) is replaced with the MGF given by ^. In ( |29] l, 
F^ (•) denotes the Lauricella Fd hypergeometric function 
of L variables, defined in 133] and 1^ Eq. (18)]. For the 
equal correlation model, the average error probability can be 
calculated using ll22l Eq. (32)], which is a special case of ( |29] l. 
Note that d29T l is here presented to demonstrate existence of 
an exact closed-form solution, whereas (|28] ) is much more 
convenient for accurate and efficient numerical integration. 
For example, the average error probability may be obtained 
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Fig. 6. Exact and approximate eiTor performance of an EGC receiver with 
correlated Nakagami-m, branches, when correlation is described by correlation 
matrices Ai and A2, rriz = 2 and /i = 



by applying the Gauss-Chebyshev quadrature rule 1341 Eq. 
(25.4.38)] over ( |28] ). In the case of the balanced diversity 
branches with equal or exponential correlation, the combina- 
tion of this quadrature rule with Tables I and II allows efficient 
and extremely accurate evaluation of the EGC performance. 

The average error probability of correlated Nakagami- 
m EGC system with BFSK modulation / non-coherent de- 
modulation is approximated by known expression of the 
equivalent MRC system 126] Eq. (16)], Pegc-bfsk ~ 
Pmrc-bfsk = \ Mr2(-\), where Mr2{-) is given by ^. 

C. Validation via Monte-Carlo simulations 

Next, we illustrate the tightness of the error performance of 
an correlated Nakagami-rn EGC system with BPSK modula- 
tion / coherent demodulation to that of the equivalent MRC 



system. The results for the actual EGC system are obtained 
by Monte-Carlo simulations, whereas those of the equivalent 
MRC system are obtained using ( |28] l. 

1) Equal and exponential branch correlation: Figs. 4 and 5 
displays the comparative error performance of the actual EGC 
and the equivalent MRC systems, for several combinations 
of (p, L, raz, ^k)- In order to accommodate unequal average 
branch powers (thus, unequal average branch SNRs), we used 
the exponentially decaying profile, modelled as 



rjfc = rjicxp(-^(fc- 1)), i<k<L, 



(30) 



where J7i is the average power of branch 1 and /i is the 
decaying exponent, with /i = denoting the case of branches 
with equal power (i.e., the balanced branches). 

2) Arbitrary branch correlation: Fig. 6 depicts the com- 
parative error performances of the EGC with same correlation 
matrices from Section II. D, Ai and A2, and the equivalent 
MRC system with respective Green's matrices Ci and C2. 
The high accuracy of our approach is maintained for arbitrary 
branch correlations. 

IV. Conclusions 

A tight closed-form approximation to the distribution of the 
sum of correlated Nakagami-m RVs was introduced for the 
case of identical and integer fading parameters. The proposed 
method approximates this distribution by using the statistics of 
the square-root of the sum of statistically independent Gamma 
RVs. Examples indicate that the new approximation is highly 
accurate over the entire range of abscissas. To demonstrate 
this more rigorously, the proposed distribution is tested against 
the computer generated data by the use of the Chi-square 
and the Kologorov-Smirnov goodness-of-fit tests. In case of 
maximal correlation, the proposed distribution becomes the 
exact distribution. 

The presented approach allowed to successfully tackle 
the famous problem of analytical performance evaluation of 
an EGC system with arbitrarily correlated and unbalanced 
Nakagami-r7i branches. The significance of the presented re- 
sults is underpinned by the existence of a large body of liter- 
ature dealing with MRC performance analysis, which permits 
highly accurate and efficient EGC performance evaluation. 

Appendix A 

Using ini Eqs. (9.212 (1)) and (7.622 (1))], one has the 
following identity 



J{m,a,p,q) = 



1 



r(m) 



- — ; to; —au ] du 






2' 2 ' l + 2a, 

(A.l) 

Using I3TI Eq. (9.212 (3)), pp. 1023] with some simple alge- 
braic manipulations, the general form ( fTsT i of the coefficient 



W{2, 1, 1) can be simplified as 



W{2,1,1)^ 



r(m. 



1/2) 
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J{m^ 
J{mz 



l,a,-l,-l) 

-l,a,-l,l) 

(A.2) 



where a = ^/(l + {N - 1)^) 



Appendix B 

The unknown joint moments in (fT4b . E[Zf-^ZiZj], 
E[Z^ZfZj], E[Z^Z,Z'j] and E[Z,nZnZiZj] can be calcu- 
lated from lfT4l Eqs. (11) and (12)]. Here we derive their 
simpler and computationally more efficient alternatives. The 
alternative to |[T4l Eq. (12)] is derived directly from the 
definition of the joint moment E[ZiZ2ZzZ4\, 



/•OO /"OO /"OO /• 

E[Z,nZnZ,Zj\ = 

Jo JO Jo Jo 



^m^n^i^i 



^fz„,z„ z, z, ( z„i , Zn , Zt , Zj ) dzm dz^ dzi dzj , (B . 1 ) 

where joint pdf of four exponentially correlated Nakagami-TO, 
RVs is expressed as lfT4l Eq. (9)] 



/. 



Z^ Z^ Zi Z- 



\ Zjyi , Zji , Zj , Zj ) 
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xI„i^_i{2m,\'4;34\ziZj) exp ( - to^(V'ii2^ 

+fe4 + V'33^:,^ + i'44z])y, (B.2) 

where ^ — [tpij] is defined by ( |23] ). Now, we integrate lfT4l 
Eq. (11)] over Zm and Zj, respectively obtaining 

Z™"+^ CXp (-TO^Vll^m) Im.,~l {2m ^Ip 12 ZmZn) 



7: {m'z1pl2Zn)"^' ^ {mztpii) ''™- + ^/2) 
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and 
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exp {-mz-ip44Z^) Im,-i {2m,ip3,4ZiZj) 



2 [m,tp34Znj (m,1p44) 



Tim^rim p L + ^-TO -TO ^34 2\ .o4^ 

We then use the series expansion of the modified Bessel 
function of first kind ||3T1 Eq. (8.445)] that allows to separate 
the integrations per variables z„ and Zj, yielding ( l22b . A 
similar procedure yields to an alternative of lfT4l Eq. (11)] 
given by (ISTT i. 
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